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ABSTRACT 

The activity levels of stars are influenced by several stellar properties, such as stellar rotation, spectral 
type and the presence of stellar companions. In analogy to binaries, planetary companions are also 
thought to be able to cause higher activity levels in their host stars, although at lower levels. Especially 
in X-rays, such influences are hard to detect because coronae of cool stars exhibit a considerable amount 
of intrinsic variability. Recently, a correlation between the mass of close-in exoplanets and their host 
star's X-ray luminosity has been detected, based on archival X-ray data from the ROSAT All-Sky Survey. 
This finding has been interpreted as evidence for Star-Planet Interactions. We show in our analysis that 
this correlation is caused by selection effects due to the flux limit of the X-ray data used and due to 
the intrinsic planet detectability of the radial velocity method, and thus does not trace possible planet- 
induced effects. We also show that the correlation is not present in a corresponding complete sample 
derived from combined XMM-Newton and ROSAT data. 

Subject headings: stars: planet-star interactions — stars: activity — stars: coronae — X-rays: stars 



1. INTRODUCTION 

The possibility of interactions between stars and their 
planets, causing an activity enhancement of the host star, 
is currently a debated issue. Possible mechanisms for such 
Star-Planet Interac tion (SPI) are tida l and magnetic in- 
teraction scenarios ()Cuntz et al.ll2000l ): also, planets trig- 
gering the release of energy, which has been built up in 
tangled cor onal loops b y normal stellar activity processes, 
is possible (|Lanzall2009t ). While chromospheric and photo- 
spheric obs ervations have yielded some hints for such in - 
teractions (jShkolnik et alj|2005l [200l lLanza et al.ll2010H . 
the analysis of p ossible coronal signatures of SPI has led to 
differing results (jKashvap et al.ll2008t Poppenhaeger et al.l 
120101: IScharfll201d: iPoppenhaeeer et al.ll2011h . 

In a recent study. IScharJ (|2010f) analyzes a stellar sam- 
ple derived from archival ROSAT X-ray data and derives 
a correlation of exoplanetary mass M p with the host star's 
X-ray luminosity for planets closer to their host star than 
0.15 AU, which is interpreted as a lower floor of possible 
stellar X-ray luminosity caused by the presence of mas- 
sive, close planets. We replot the data from that sample 
in Fig. [TJ it indeed shows an amazing correlation of plane- 
tary mass and stellar X-ray luminosity. This would be an 
extremely interesting finding if the correlation was really 
caused by SPI. We therefore conduct an in-depth analysis 
of possible factors abl e to cause su ch a correlation, using 
both the sample from Scharf (2010) as well as a complete 
sample of planet- hosting stars within 30 pc d istance, which 
we presented in iPoppenhaeger et al.l ((2010) . We specif- 
ically investigate possible sample selection effects, deter- 
mine suitable variables which should be tested for correla- 
tions with each other, and compare the results derived for 
the two samples mentioned above. 

2. SAMPLE PROPERTIES 

IScharJ (|2010f) selected a stellar sample consisting exclu- 
sively of X-ray detections (and upper limits) derived from 
ROSAT All-Sky Survey (RASS) data. This sample has two 



shortcomings when compared to the wealth of data avail- 
able from t oday's X-ray missions XMM-Newton and Chan- 
dra: First, ISchard (l2010D uses no pointed observations, ar- 
guing that in RASS a given star was scanned several times, 
yielding a time-averaged X-ray luminosity. Howev er, this 
is only true for stars close to the ecliptic poles, as IScharj 
(|2010[) points out correctly. Stars near the ecliptic equator 
are nominally scanned several times only over a two-day 
period, and most of the given stars with close-in planets 
are located at latitudes of ±45° or lower. Additionally, 
the orbital periods of planets in the given close-in sub- 
sample range from 2.2 — 8.7 d, yielding phase coverage of 
substantially less than one orbit for almost all cases. Thus 
there is no specific advantage of the RASS data compared 
to XMM-Newton or Chandra data, especially given the 
higher sens itivity and spatial accuracy of these two X -ray 
telescopes (jJansen et al.ll200TI : IWeisskopf et al.ll2000[) and 
the short total observation time of RASS sources which 
is usually well below 1 ks, compared to typical exposure 
times in pointings which are in the 10 — 100 ks domain, 
depending on the target. 

Second, the two stars in the close-in sample with 
the highest X-ray luminosities, i.e., HD 41004B and 
HD 162020, are not comparable to the rest of the sample. 
HD 41004AB is a binary system consisting of a K star 
and an M star in a close orbit with a projected distance 
of 0.5", the M star being the component with the sub- 
stellar companion which is at the boundary of planet and 
brown dwarf. The RASS data does not allow to determine 
the X-ray luminosity of each of the two stars individually 
because of ROSATs rather broad FWHM; however, M 
stars us ually have lower averag e X-ray luminosities than 
K stars (|Schmitt fc Liefkd [20041 ) . so that the complete X- 
ray luminosity of the system should not be attribut ed sole- 
ley to the M star with the planetary companion. IScharJ 
(|2010f) notes that unresolved binarity contributes to the 
errors in his sample, however in this case the binary na- 
ture of HD 41004 is known and can be accounted for. The 
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Fig. 1. — Left: Lx vs. planetary mass for stars with close-in planets, data from Scharf (2010). The two incomparable data points form 
the original sample are shown as red solid squares. For comparison, X-ray data for 51 Peg is inserted as blue open triangles. Right: Activity 
indicator Lx/Lb i vs. planetary mass for stars with close-in planets, data from Scharf (2010). 



other sta r, HD 162020, mi ght be a young star as is dis- 
cussed in lUdrv et al.l (|2002l ); however, the age determina- 
tion for this star is not entirely clear. If it is a pre-main 
sequence star, it would not be comparable to the rest of the 
sample, since such stars are known to have much higher 
X-ray luminosities than their main-sequence counterparts 
(jFeigelson fc Montmerlej|1999l) . 

3. METHODOLOGICAL CAVEATS 

When testing for correlations of stellar activity with 
other quantities, an appropriate proxy for stellar activ- 
ity has to be chosen. Choosing X-ray luminosity as this 
proxy has two disadvantages: First, there is a correla- 
tion between stellar radii and stellar X-ray lum inosities, 
at least for stars with outer convection zones (|Schmittl 
Il997t iRaassen fc Kaastral l2006t). This is why one usually 
normalizes the X-ray luminosity with the bolometric lumi- 
nosity of the star to make the activity levels of stars with 
different spectral types comparable. As a rule of thumb, 
stars with log Lx jLboi < —5 are dubbed "inactive", while 
stars with log LxjL^oi > —4 are dubbed "active", regard- 
less of spectral type. So, an analysis of X-ray activity for 
a sample with a variety of stellar masses should also test 
for correlations of planetary parameters with Lx / L^oi to 
check whether some stellar mass distribution causes a fake 
signal in correlations with Lx- 

Second, in flux-limited survey data there is usually a cor- 
relation between detected luminosities (and upper limits) 
and the distance d of the targets. This is a consequence of 
the approximately constant exposure time in a survey and 
the strength of the source signal of targets which scales 
as 1/d 2 . So, when dealing with survey data, one should 
carefully check for dependencies on distance. 

4. RESULTS 

The considerations above lead to two compleme ntary 
analyses of the results presented in IScharJ (|2010f ): re- 
analyzin g the origin al ROSAT sample as given in Ta- 
ble 1 in lScharj (|2010l ) for dependencies on distance d and 
Lx/Lb i, and testing for such dependencies in the much 



larger complete sample us ed in our previous SPI study 
(jPoppenhaeger et al.ll2010T ). 

4.1. The RASS data revisited 

In Fig. [TJ we show the X-ray luminosities of planet- 
hosting stars with a p < 0.15 AU from the R OSA T sample 
as a function of the innermost planet's mass. The data 
points for HD 41004B and HD 162020 are plotted in red 
for comparison. We would like emphasize the special case 
of 51 Peg, a star not included in the original sample: this 
star has been observed in several pointed X-ray observa- 
tions with ROSAT, XMM- Newton and Chandra, covering 
different phases of the planetary orbit; it remained unde- 
tected i n the RASS data and was therefore included by 
IScharj (pOlOh only as an u pper limit. A detailed an aly- 
sis of these data has shown (|Poppenha,ger et al~| [2009) the 
star's X-ray flux to be constant and at a very low level over 
16 years, indicating that the star might be in a Maunder 
minimum state despite its close-in heavy planet. This sys- 
tem is a sig nificant outli er of the Lx vs. M p i relation 
presented in lScharJ (|2010D (see Fig. [U), although it fulfills 
the criterion of having a determined phase-averaged X-ray 
luminosity. 

To test for dependencies on the distance d of the 
stars, we condu ct a Principal Component Analysis (PCA) 
( Pca rsonlll901f ) on the three lo garithmic var iables log Lx , 
log M p and logd from the full Scharf (2010) sample. We 
use the two eigenvectors with the largest eigenvalues as the 
feature vectors; the third eigenvalue is very small by com- 
parison (0.05 vs. 0.95 and 0.14), meaning that one loses 
only very little data variance in this analysis. After repro- 
jecting the data, all three variables show a strong linear 
trend with respect to each other (see Fig. [2]). This is a clear 
indication that the stellar distance is a crucial parameter 
in this sample which cannot be ignored. It is important 
to note that in an unbiased sample, Lx and M p must not 
depend on stellar distance; if they do as in this sample, a 
selection effect is present. This provides an explanation for 
the apparent dependence of Lx on M p : the detection limit 
of Lx increases with increasing distance. The detectability 
of planets is somewhat intricate, and we investiagte depen- 
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Fig. 2. — Principal Component Analysis results for the Scharf (2010) sample, stars with close-in planets given as filled symbols. All three 
parameters (d, Mp, Lx) show linear trends with respect to each other, indicating that the stellar distance has a crucial influence in this 
sample. 



dencies in detail in section [5] In short, a dependency of 
planetary mass on stellar distance is present. So, at larger 
distances, the radial velocity method favours the detection 
of heavier planets, and low X-ray luminosities cannot be 
detected any more in the R OSA T All- Sky Survey, which 
yields the observed trend of Lx with M p without having 
to invoke effects from supposed Star-Planet Interactions. 

Also without performing a PCA on this sample, the de- 
pendencies of Lx and M p on d are revealed by rank corre- 
lation tests. We calculate Spearman's p, a rank correlation 
coefficient which displays a perfect correlation by a value 
of 1, perfect a nticorre l ation by —1 and no correlation by 
0. For the full IScharJ (12010D sample, we find strong cor- 
relations of both Lx and M p with d, indicated by p val- 
ues of 0.49/0.54 respectively, translating to probabilities of 
0.6/0.2% that such a correlation can be reached by pure 
chance. This correlation analysis yields the same result as 
the PCA; the stellar distance is the crucial parameter in 
this sample which causes the Lx /M p correlation. 

This is also reflected in the behavior of Lx/Lf, i, where 
there is no significant correlation with planetary mass for 
stars with close-in planets, see Fig. Q] We also checked 
this visual result with a Spearman's p test while exclud- 



ing the data from the two incomparable stars. This yields 
p = 0.05, i.e. a very weak positive correlation; the proba- 
bility that such a p value is reached by chance is 87%. This 
is not surprising: if the trend in Lx is a distance selection 
effect and not related to the stellar activity level, then 
the quantity Lx / L^d , which measures the intrinsic stellar 
activity level, should be independent from the planetary 
mass. 

4.2. The correlation as seen with XMM-Newton 

In our further an a lysis, we use the data presented in 
iPoppenhaeger et al.l (|2010l) . which consists of all known 
planet-hosting stars within a distance of 30 pc from 
the Sun, with X-ray properties derived from XMM- 
Newton and ROSAT dat a. As discussed in detail in 
IPoppenhaeger et ail (|2010[) . we exclude unresolved binary 
stars and early stars without outer convection zones from 
our analysis. The errors given are Poissonian errors plus 
an additional uncertainty of 30% on the X-ray luminosity 
to account for short-time variations, since a large part of 
our sample consists of pointed XMM-Newton obervations. 
We use the sa me samp l e sele ction criterion on these data 
as was used in IScharfl ()2010[ ) (planets at a < 0.15 AU). 
We show the relation between Lx and M p i in Fig. [3J left 
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Fig. 3. — Left: Lx vs. planetary mass for stars with close-in planets (a < 0.15 AU), data from Poppenhaeger ct al. (2010). The higher 
sensitivity of XMM-Newton yields many a dditional X-ray detections in the lower right corner of the diagram, compared to Fig.[T] Stars which 
are also present in the sample from[Scharf (20 lT|) are plotted as gre en filled symbols. Right: Activity indicator Lxl^bol vs. planetary mass 
for stars with close-in planets, data from [Poppenhaeger ct al. (20ljj). No significant correlation is present. 
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Fig. 4. — Same as Fig. [3] but for stars with planets at semimajor axes > 0.15 AU. 



panel; d ata from stars which are also present in the sam- 
ple from Scharf (2010) are plotted as green filled symbols. 
These stars lie close to a straight line, similar to Fig. [TJ 
although the data was collected in single pointings and 
not averaged over larger portions of the planetary orbit. 
This shows that the averaging process is not crucial for 
this kind of analysis; the Lx values derived from XMM- 
Newton pointings are very similar to the ones from the 
RASS data. The main difference to Fig. []] is that there 
are many additional X-ray detections in the lower right 
corner of the diagram. This is contrary to the assumption 
that massive, close-in planets cause a lower floor for the 
X-ray luminosity of their host stars. 

Additionally, in this sample there is no significant cor- 
relation in the relation between the X-ray activity indica- 
tor Lx/Lboi vs. planetary mass (see Fig. [3j right panel); 
testing for rank correlation yields p — 0.003, i.e. practi- 
cally no correlation at all. This is also true for stars with 
far-out planets, for which no SPI-related effects are ex- 
pected (Fig. 2]) . The only significant correlation present in 



the whole sample is one between X-ray luminosity and the 
product of planetary mass and inverse semimajor axis. For 
the intrinsic X-ray activity measured by Lx/L^pi no such 
correl ation is present. As discussed in lPoppenhaeger et al.l 
(|2010l ). the Lx correlation is equally strong in a subsam- 
ple of stars with small, far-ou t planets as well as in stars 
with heavy, close-in planets. iPoppenhaeger et ail (j2010) 
conclude that the correlation is caused by selection effects 
from planet detection; if it was caused by SPI, it should 
be strong in the second subsample and weak in the first 
subsample. 

For the sake of completeness, we also conducted a PC A 
for this sample. Here we find for the reprojected data that 
there is a strong linear trend of M p with d, but no appar- 
ent trends of Lx with d or M p . This is due to the fact that 
we also use data from pointed observations in our sample, 
where observations of more distant targets usually have 
longer exposure times. This prevents the correlation o f 
Lx with d which is present in the sample of lScharj (|2010f ). 
and therefore also no correlation between Lx and M p is 
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present here. 

4.3. The nature of the correlation between distance and 
planetary mass 

In the above section we used a volume-limited sample 
of planet-hosting stars in which most of the stars have 
been detected in X-rays. In that sense, the sample is 
complete with regard to X-ray flux. However, the sam- 
ple is most probably not complete regarding the detection 
of planets. Therefore, selection effects induced by plane- 
tary detectiojijujrthojfe carefully (see 
also lOToole et al.l (|2009l ) : iHartmanl (|2010D ). The portion 
of planets which are detected by transits is quickly grow- 
ing since the start of the CoRoT and Kepler observations, 
but for stars in the solar neighborhood, the dominant de- 
tection mechanism still is the radial velocity method. We 
now provide an investigation of such selection effects in the 
basic planetary and stellar parameters that are present the 
sample used above. 

The RV-detectability of a planet depends on several 
properties: first, on the brightness of the star itself and 
therefore the quality of the signal in which one searches 
for RV variations; and second, on the stellar mass M*, the 
planetary mass M p , and the planetary period P. Other 
influences such as eccentricity of the planetary orbit are 
ignored here. Specifically, the RV semi-amplitude is pro- 
portional to P- 1 / 3 x M p x MT 2/3 . Thus it should be 
easier to detect low-mass planets around low-mass stars 
for a given (fixed) sensitivity. 

The first half of our hypothesis is that the RV detectabil- 
ity decreases as the stellar distance d increases, since the 
apparent brightness of the star decreases. This means that 
at larger distances, planets should be found around intrin- 
sically brighter stars. On the main sequence, this implies 
earlier and therefore more massive sta rs. We test this with 
Spear man's p for the sample from iPoppenhaeger et al.l 
(2010); we use the complete sample now, not only the 
subsample with close-in planets, as we aim to analyze se- 
lection effects for all detected planets. We find that indeed 
the stellar mass is strongly positively correlated with dis- 
tance (p = 0.36, Pfalse < 1%). 

The second half of our hypothesis is that because of 
the dependencies of the RV amplitude on stellar and plan- 
etary parameters, mostly massive planets in close orbits 
will be found around massive stars, as the RV amplitude 
will be too small to be detected otherwise. Accordingly, 
wc find that stellar mass and planetary mass are positively 
correlated (p = 0.46, Pf a ise < 0.1%); the correlation be- 
tween stellar mass and the quantity P" 1 / 3 x M p is even 
stronger with p — 0.53). We also show the correlation of 
M* and P -1 / 3 x M p in Fig. [5] However, for a linear graph- 
ical representation, one has to choose specific functions 
(for example, logarithms) of the parameters to be plotted. 



This may not be the real relation of the two quantities; 
furthermore, we have already shown that there is an addi- 
tional dependence of the RV detectability on stellar mass 
and distance which may also influence other correlations. 
Therefore, this plot in which a linear trend is indicated 
does not directly represent the true extent of the corre- 
lation as demonstrated by the rank correlation analysis. 
Indeed, it is advantageous to use rank correlations here 
since they do not rely on any specific parametrization of 
the quantities in question. 

The above analysis confirms our detectability consider- 
ations. At larger distances, planets are detected around 
brighter, more massive stars, and these planets are more 
massive and in closer orbits in order to produce a strong 
and therefore detectable RV signal. 

In our PCA calculations in section 14.11 (Fig. the 
planets in wider orbits, depicted by open symbols, show 
a larger scatter around the various correlations than the 
close-in planets. When we analyze only the systems with 
wide-orbit planets (a > 0.15 AU) for rank correlations, we 
find that the correlations between planetary mass and stel- 
lar mass as well as between stellar mass and P- 1 / 3 x M p 
are of comparable strength as in the full sample. However, 
the correlation between distance and stellar mass is weaker 
(Pfaise ~ 10%). One could speculate that in the outer- 
planets sample the RV signal is weaker due to the larger 
distance between star and planet, and this might favour 
detections around brighter stars, yielding less "leverage" 
for the correlation (the mean stellar mass is actually larger 
for the outer-planet-sample). However, we do not have a 
compelling and easy explanation for this. 

In summary, we deduce that the correlation of Lx and 
Mp is a combined selection effect of X-ray flux limits in 
the ROSAT All-Sky Survey and planet detectability by 
the RV method. In the sample that is not X-ray flux lim- 
ited, the correlation of planetary mass and X-ray activity 
is consequently not present, as demonstrated in Fig. [3] 

5. CONCLUSIONS 

We conclude that there is no detectable influence of 
planets on their host stars, which might cause a lower floor 
for X-ray activity of these stars. Rather, possible planet- 
star interactions seem to induce only small effects on the 
host stars, which will however provide valuable informa- 
tion on stellar and planetary magnetic fields if measured 
in X-rays. 

We thank an anonymous referee for his/her thoughtful 
comments and the detailed discussion of the analysis pre- 
sented in our paper. K. P. acknowledges financial support 
from DLR grant 50OR0703. 
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